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ABSTRACT 

Galaxy groups are quite underluminous in X-rays compared to clusters, so the intracluster medium 
has to be considerably underdense in the former. We consider this to be due to substantial energy fed 
back into the ICM when the baryons in the member galaxies condense into stars ending up in SNe, 
or accrete on to central supermassive black holes energizing AGNs. We compute the outflow and the 
blowout effects driven by the AGNs and the resulting, steep luminosity-temperature correlation Lx — T. 
We compare this with the SN contribution and with the X-ray data; the latter require the AGN energy 
to be coupled to the surrounding ICM at fractional levels around 5 10~^. We link the Lx — T behavior 
with the parallel effects of the AGN feedback on the gas in the host galaxy; we find that these yield a 
correlation steep up to M, oc between the galactic velocity dispersions and the central BH masses. 

Subject headings: galaxies: clusters: general - quasars: general - shock waves - X-rays: galaxies: clusters 



1. INTRODUCTION 

Groups and clusters of galaxies shine in X-rays due to 
the thermal bremsstrahlung radiation emitted by the hot 
intracluster medium (ICM) they contain. 

But the poorer groups are found to be progressively un- 
derluminous, so as to lie substantially below the simple 
scaling Lx oc T^. For the latter to hold, in the luminosity 
Lx oc 7^1/2 ICM number density n would have to 

be proportional to the gravitationally dominant DM mass 
density p. This constraint adds to the temperature T be- 
ing close the virial value Ty and the size R scaling as the 
virial radius Ry oc Ty^^ 

In fact, the observed Lx — T correlation has a shape 
more like Lx oc for richness 1 clusters (Kaiser 1991), 
and in moving toward poor groups it bends further down 
to Lx oc Ty or steeper (Ponman, Cannon, & Navarro 1999; 
see also § 3). So in groups the ICM is considerably under- 
dense relative to the cluster values n w 10~^ p/TTip- 

Correspondingly, the central ICM entropy S/k = 
In kT /n?/^ deviates upwards from the simple scaling e'^ oc 
Ty^ to attain in poor groups the 'floor' value kTjr?!^ « 
140 keV cm^ (Lloyd-Davies, Ponman & Cannon 2000). 
This requires a density deficit associated with increased 
or constant T/Ty. Such a non-adiabatic behavior may be 
traced back to energy losses from, or additions to the ICM. 

Here we focus on the 'heating' that arises as the baryons 
in the member galaxies condense into stars followed by 
SNe, or accrete onto a supermassive black hole (BH) kin- 
dling an AGN. Two issues stand in the way: (i) any extra 
energy from sources has to compete with the huge thermal 
value E lO^^ {kTy/\eNf/'^ erg for the ICM in equilib- 
rium; (ii) the extra energy that can be coupled to it is still 
poorly known. We discuss these issues, and derive two 
observables that bound or probe at group and at galactic 
scales the amount of energy coupled. 



2. FEEDBACK FROM SNE 

Obvious first candidates for energy discharges into the 
ICM are the SN explosions following star formation in the 
member galaxies of groups and clusters. Prompt, type II 
SNe canonically release 10^^ ergs; these are effectively cou- 
pled to the gas when SN remnants propagate cooperatively 
over galactic scales to drive galactic winds (Ostriker & Mc- 
Kec 1988; Wang et al. 2001; Pettini et al. 2001; Heckman 
2002). With a coupling around 1/2, the energy input (in- 
cluding winds from hot stars) comes to < 3 lO''^ erg 
per solar mass condensed into stars. In a fiducial group 
with kTy w 1 keV, virial mass M„ w 5 10^^ M© and stel- 
lar mass around 3 10^^ Mq, this would raise by kAT « 0.3 
keV the temperature of the entire ICM. The outcome looks 
like a modest AE/E = AT/Ty « 1/3. 

Actually, SNe make optimal use of their energy in that 
they produce hierarchical preheating of the ICM, while 
a group and its ICM are built up hierarchically through 
merging events with a range of partners. About half the 
final DM mass My is contributed to the main progenitor 
by smaller partners with masses within the window My/S 
to My/20, corresponding to temperatures from 0.6 down 
to 0.15 (Lacey & Cole 1993; Menci & Cavaliere 2000). 

Smaller lumps in the window have shallower gravita- 
tional wells and produce more star-related energy on scales 
closer to the dynamical time td] so they are more effective 
in heating/ejecting their gas share. During each subse- 
quent step of the hierarchy such gas preheated externally 
will less easily flow into the main progenitor's well. Thus 
the effects propagate up the hierarchy, and lower ICM den- 
sities are induced in all structures up to poor clusters. 

Two main density suppression factors arise in moving 
from clusters to groups. These are best discussed on re- 
ferring to ICM in hydrostatic equilibrium within the DM 
potential well Acj) (Cavaliere & Fusco-Femiano 1976; Jones 
& Forman 1984). In the corresponding ICM density run 
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n{r) = n2 exp [/3 A(j>{r)] the energy injection resets the 
values of (3 = Ty/T and of n2- 

One factor is related to the outflow and is given by (3; 
the latter is lowered by about 0.6 from rich clusters toward 
poor groups where stellar preheating provides a contribu- 
tion to T comparable to T^. A second factor is the differ- 
ential containment expressed by the boundary value n2. If 
this is set by jump conditions across the accretion shocks 
at r ~ (see Takizawa & Mineshige 1998; Gheller, Pan- 
tano & Moscardini 1998), the density is further suppressed 
from clusters to groups by another factor approaching 1/2 
on average, see Menci & Cavaliere (2000). 

These authors model the process on grafting the ICM 
equilibrium onto the semi-analytic treatment (SAM) of 
star and galaxy formation. This is based on the merg- 
ing histories of the DM, and includes star formation and 
gas heating/ejection by SNe in terms of simple recipes. 
The latter imply heating to dominate at scales of bright 
galaxies and larger; only in small galaxies the gas fractions 
blown out exceed 10~^ (sec Madau, Ferrara & Rees 2001). 

With SN feedback, the SAMs provide detailed flts to 
the stellar observables. But the agreement with the avail- 
able Lx — T data for poor groups is marginal if both the 
Navarro, Frenk & White (1997) potential and the stan- 
dard flat A-cosmology are adopted. This is illustrated in 
next § 3, see also Borgani et al. (2001). 

Note that the SAMs include the amount of cooling suit- 
able for baryon condensation to stars, which helps in low- 
ering the ICM density. A leading role of cooling (combined 
with suitable feedback) to remove the low entropy gas has 
been discussed by Bryan (2000), Muanwong ct al. (2002), 
and Voit et al. (2002). We will concentrate on a maximal 
but realistic feedback process. 

3. FEEDBACK FROM AGNS 

The natural sources of strong feedback are the AGNs, 
energized by accretion of cool baryons onto supcrmas- 
sive BHs in galactic cores (see Wu, Fabian & Nulsen 
2000; Bower et al. 2001). The outputs are large, of or- 
der 2 10*^2 (M,/10^ Mq) ergs for an accreted mass M,, 
with the standard mass-energy conversion efficiency of or- 
der 10^^. If a fraction / is coupled to the surrounding 
medium, the energy actually injected comes to Ai? « 
/ 10^" erg per solar mass condensed into stars; we have 
used 2 10"^ for the ratio of the BH mass to that of the 
current host bulges (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000; see also Fabian & Iwasawa 1999). 

Compared with SNe, AGNs potentially provide a larger 
energy output in shorter times, close to td of the host 
galaxies. However, / is more uncertain than the analo- 
gous quantity for SNe, though is expected to be lower. 

The 10% radio-loud AGNs directly produce large kinetic 
energies in the form of jets, but to now the observations 
indicate limited impact on the gas surrounding a number 
of active sources (see McNamara et al. 2001; Terashima 
& Wilson 2001; Young, Wilson & MundeU 2002). In the 
90% radio-quiet AGNs a smaller coupling / « 10~^ is ex- 
pected in view of their ffat spectra and of the low photon 
momenta. Values up to / ~ 10~^ are only conceivable in 
systems where the photons are heavily scattered / absorbed 
within the gravitational reach of the BH, and escape in 
hard X-rays if at all (Fabian, Wilman & Crawford 2002). 



Whence the interest in probing / from overall effects on 
the ICM. During the AGN activity the gas initially con- 
tained in a large galaxy or a group will be heated up and 
partially blown out. We will treat blowout and outflow 
driven internally in poor groups with kT^ ^ 1 kcV where 
AE/E ~ 1. These effects are maximized in spherical sym- 
metry, not unflt for the radio-quiet AGNs we consider. 

Within the structure's dynamical time tj, we describe 
the transient regime as a blastwave sweeping through 
the surrounding gas (see Platania et al. 2002); when 
AE/E < 1 holds the leading shock at r = i?s is not nec- 
essarily strong, and DM gravity is important. 
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Fig. 1. — Radial distributions of density and temperature (nor- 
malized to their postshock values) in a blastwave driven by a flaring 
AGN. Dot-dasfied and dotted lines refer to the temperature behind a 
strong {AE/E = 1.8) or a weak shock (AE/E = 0.3), respectively; 
solid and dashed lines refer to the density. On approaching the pis- 
ton (r = 0.66 -Rs or r = 0.45 -Rs for the adopted values of AE/E) 
the density diverges weakly while the mass inside r vanishes. 

From the relevant hydrodynamical equations we have 
derived (and illustrate in fig. 1) a family of self-similar 
solutions of the Sedov (1959) type that include the DM 
gravity, a steep initial density gradient, and centrally in- 
jected energy AE(t) growing over times of order t^. The 
radiative cooling is slow on mass- average in our groups. 

Our fiducial case will have AE{t) oc t injected into an 
initial configuration with n(r) oc r~^, i.e., isothermal ICM 
in hydrostatic equilibrium (/3 « 1) in the potential pro- 
vided by DM density p{r) oc r~^; we denote by E{Rs) the 
modulus of the total initial energy within the shock radius 
Rs- Then the leading shock moves outward at a costant 
speed Vs, only moderately supersonic when AE/E < 1. 

Self-similarity implies AE{t) / E{Rs) to be independent 
of time and position, as is especially simple to see in our 
fiducial case where E{Rs) oc Rg oc t oc AE{t). For two 
values of AE/E we show in fig. 1 the density and tem- 
perature runs. The flow begins at a 'piston', the inner 
contact surface where the density diverges weakly while 
the gas mass within r and the temperature T(r) vanish. 

In fact, the perturbed gas is confined to a shell with 
inner (piston) radius XRg and outer (shock) radius Rg. 
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Self-similarity implies the thickness ARg/Rg ~ 1 — A of 
such a shell to depend only on AE/E; for strong shocks 
driven by AE/E > 1 we find A 0.84, while for a weak 
shock corresponding to AE/E = 0.3 we find A = 0.45. 

These considerations lead us to represent our solutions 
in terms of the simple 'shell approximation', known to pro- 
vide results reliable to better than 15%, see Cavaliere & 
Messina (1976); Ostriker & McKee (1988). In this approx- 
imation the energy balance reads 



AE + E = 



1 2 3 GMm 



(1) 



and directly shows the relevance of AE/E. Here M is the 
DM mass within Rs\ also, V = AttRI (1 - A^)/3 is the vol- 
ume of the shell; m and p are the associated gas mass and 
mean pressure; finally, V2 oc Vg is the postshock velocity 
given by the Rankine-Hugoniot conditions. 
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Fig. 2. — The Lx — T correlation; bolometric Lx including 
standard line emissions. Thin dotted line: gravitational scaling 
L oc T2. The shaded one-sigma strip results from a SAM including 
the stochastic merging histories of the DM and the SN feedbax;k, 
see § 2. Thick lines: our results with feedback from AGNs accret- 
ing as specified at the end of § 3, and with / = 3 10^"^ (solid) and 
/ = 10~i (dashed). Data: Markevitch (1998, circles), Arnaud & 
Evrard (1999, squares), Helsdon & Ponman (2000, triangles). 



Self-similarity requires all terms in eq. (1) to scale like 
Rs] the coefficients depend only on AE/E, and are eas- 
ily derived following the pattern indicated by the above 
authors. We find the ratio of the kinetic to the thermal 
energy (i.e., the 1*^' to the 2"'^ term on the r.h.s. of eq. 1) 
to range from 5 lO^'^ up to 2 when AE/E increases from 
0.3 to values larger than 1. Analogously, one may derive 
the dependence of Vg on AE/E, and an analytic approxi- 
mation to the mass distribution within the shell. 

After the passage of the blastwave, and before a ma- 
jor merging event reshuffies the DM mass substantially, 
the gas recovers hydrostatic equilibrium and again n(r) = 
exp [(3 A0] holds; but now the governing parameters are 
those given in Table 1. The value oi (3 = T^/T (related to 
outflow) is reset using the mass-averaged temperature T. 
The new ICM mass m — Am (left behind by the blowout) 
is taken to be that still residing at t = td between the 



piston and Ry; thus the boundary condition n2 is reset by 
requiring consistency with the volume integral of n(r). 

We then compute Lx oc T^/^ J dr n'^{r) and plot it as 
a function of temperature in fig. 2; here we approximate T 
with Ty, since these differ only modestly as seen from the 
values of (3 in Table 1. Our results are given for two val- 
ues of the energy coupled; these bracket AE = 2 10®*^ erg 
corresponding, e.g., to / w 510^^ and to M, = li]'^ Mq 
for the largest BH (or sum of BHs) formed within td in 
groups with membership around 10 bright galaxies. 

The variable Ty may be related to the quantity AE/E 
that governs the blastwave on using E oc {kTy)^/"^, see § 1. 
Assuming AE oc M, (see the beginning of this Sect.), the 
correspondence is simply given by 

E I 10-2 y I 109 M^y \keV J 



(2) 



In fig. 2 we have actually implemented the second approx- 
imation AE oc M. [1 -h (fcT„/keV)6/^]/2. 

This accounts for the cosmological evolution of the 
AGNs that occurs at redshifts z < 2.5, when groups and 
clusters with increasing kTy oc My^^ (1 -I- z) are formed 
by the standard hierarchical clustering. In the critical 
universe the above (normalized) correction stems from: 
(a) the diminishing fraction of member galaxies activated 
within the increasing td, that scales simply as My, i.e., 
approximately as (1 + z)~^; (b) all outputs weakening as 
(1 -|- z)^ (see Cavaliere & Vittorini 2002). In the standard 
fiat A-cosmology the result provides a close upper limit. 

In fig. 2 we also recall the contribution from SNe, and 
report the available X-ray data. 

Table 1 

Parameters of the recovered equilibrium 



AE/E 


/3 


1 — Am/m 


0.3 


0.94 


0.92 


1 


0.86 


0.58 


1.8 


0.8 


0.16 



* The approximation Am/m ~ 0.5 AE/E holds to better 

than 10% for AE/E < 1.4. 

4. DISCUSSION AND CONCLUSIONS 

Two key features are apparent in fig. 2 and are spelled 
out in Table 1. First, into the cluster range the devia- 
tions from the gravitational scaling vanish because both 
1 — Am/m and (3 = Ty/T saturate to 1. Second, moving 
into the group range the luminosity is non-linearly sup- 
pressed as Lx oc n| oc (1 — Am/m)"^, due to the increasing 
contribution from the blowout as AE/E raises toward 1. 

The current X-ray data in groups are seen to require 
values around / « 5 10-^ in our blastwave model. With 
these, the feedback from AGNs dominates over SNe in 
poor groups, causing stronger suppression of Lx and fur- 
ther bending of the Lx — T relation. Variance of / from 
3 10^2 to 10"^ produces a widening strip, but one still con- 
sistent with the current data and their scatter. This also 
covers the effects of moderate non-spherical deviations. 

As our Lx — T relation intrinsically steepens towards 
poor groups, we can check it in the adjoining galactic range 
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where cooling still does not dominate. In terms of the ve- 
locity dispersion a = {kTy / 0.6 nip)^^^ from the virial re- 
lation, we find a steepening correlation Lx oc cr"; for the 
upper values of / the slope is n « 8.5 10 in large galaxies 
with (T = 300 km s^^, which accords with the detections 
and the fit by Mahdavi & Geller (2001). 

Our result is due to a blastwave with AE/E w 1.2 caus- 
ing Am/m « 0.6. Down to what scales can we extend the 
increasing trend of AE/E? We argue AE/E on average 
will not excecid 1 by much, nor will Am/m attain 1. 

First, let us impose the limiting constraint AE/E w 1 
to eq. (2), and find the accreted masses 

M. »2109Mq f-^^ ( j\ . (3) 

" VlO-V V300kms-V 
Such values, consistent with those adopted in our compu- 
tations, for / in the range 3 — 5 10^^ gratifyingly agree 
with the masses of dark objects detected at the center of 
many galaxies. Also the trend accords with the similar 
correlation pointed out by Ferrarese & Merritt (2001) and 
Gebhardt et al. (2001). Note that on using the second 
approximation to Af, discussed after eq. (2), our correla- 
tion is somewhat flatter than M, cx it has the slope 
4.3 around a w 300 km s~-^, and the prefactor 3 10^ Mq. 

Next we discuss how our limiting value AE/E « 1 arises 
in galaxies from BH accretion regulated by the AGN itself 
(sec also Silk & Rces 1998). On the one hand, sustaining 
AE/E to about 1 requires sufficient cold gas made avail- 
able for inflow. The requirement is met by gravitational 
torques exerted in the host by companion galaxies within 
small groups, see Cavaliere & Vittorini (2002) . During en- 
counter or flyby times of order td such interactions desta- 
bilize fractional gas masses of order 10~^, while the values 
needed to satisfy eq. (3) in the host galaxies are only of 
order a'^/frjc^ w 2 10"^ (cr/SOOkms-i)^. 



On the other hand, AE/E will be limited if the accre- 
tion itself can be affected on the time scale td by the AGN 
feedback; if so, the AGN will fade out. Declining luminosi- 
ties are included in our self-similar blastwave family under 
the form L cx t^i^-^)/^ if the initial density gradient fol- 
lows a steeper law n oc r~'^ with w > 2. Increasing w up 
to 2.5 corresponds to L{t) going from constant to a spike; 
up to w « 2.4 the non-linear behavior of Lx — T around 
Am/m 1/2 is generic. 

But when lo approaches 2.5 the time scales effective 
for L{t) become quite shorter than td; the correspond- 
ing blastwaves cause larger values of Am/m at a given 
AE/E. This behavior is indicative of runaway conditions 
prevailing when a galaxy happens to grow a large BH in 
short times. Then most galactic gas is blown away beyond 
Ry, so the star formation activity is suppressed already 
at 2; w 2 (see also Granato et al. 2001). Such may have 
been the case for some of the recently discovered EROs 
(see Cimatti et al. 2002; Alexander et al. 2002). 

To sum up, we find that the AGN feedback sharply 
steepens the Lx — T correlation in the poor group range, 
and links its shape to that of the galactic M, — a corre- 
lation. This is because in moving from clusters to groups 
the energy AE injected by AGNs within td grows rela- 
tive to the unperturbed E, and overwhelms SNe. But on 
entering the galactic range, AE/E approaches unity and 
constrains the accretion itself. These correlations provide 
two linked but obscrvationally independent probes of the 
hidden parameter /. On the basis of our blastwave model, 
the existing X-ray data concur with the optical ones to 
indicate values around / « 5 10~^. 

We thank G. De Zotti for stimulating discussions and 
our referee for constructive comments. 
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